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ABSTRACT 



Turbulence data obtained over natural water waves were analyzed using 
joint probability distribution and conditional means methods. These data 
represented conditions when the waves were decaying and when the waves 
were building. In both cases, significant wave-induced fluctuations were 
identified in the airflow. All features of the velocity fluctuations 
were examined for two levels above mean water level. In the case of a 
decaying wave field, decelerations in the airflow can be associated v/ith 
an assumed propagating pressure maximum over the crest of the wave. Other 
than this deceleration, the airflow appears to reflect simple streamline 
bending over the mobile irregular wave surface. In the case of a building 
wave field, velocity fluctuations appear to agree with those predicted by 
linear wave generating theories. 
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I. INTRODUCTION 



This study examines turbulence data obtained over natural water waves 
in order to determine the influence of the waves on the overlying air. 

The data examined are simultaneous measurements of velocity and temperature 
fluctuations at two fixed levels above mean water level. 

The study extends a previous analysis on data for which significant 
evidence of wave-induced motion in the airflow was obtained. Joint proba- 
bility density and conditional mean analysis are the statistical procedures 
used to extend the previous interpretations. 

The lack of good descriptions of the interaction between the ocean 
surface and the adjacent air layer is a major obstacle preventing solution 
of a wide class of applied geophysical problems such as weather prediction 
and general circulation studies. A suitable description requires the 
application, to observational data, of specific statistical methods which 
could reveal properties of wave-induced organized motion. 

Current problems requiring additional knowledge on the dynamics of the 
surface layer are the determination of: 

1. kinetic energy dissipated by boundary layer turbulence; 

2. thermal energy gained by turbulent transport of sensible and 
latent heat; 

3. energy transformed by vertical velocities due to the horizontal 
distribution (convergence) of surface friction; and 

4. surface winds used in semi -empirical wave forecasting models. 

Although numerous observational investigations have been conducted to 

describe the structure of the near surface layer, the majority of these 
measurements have been over land. Therefore, parameterization of the 
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boundary layer over the sea has been based on overland results. It is 
clear, however, that a mobile irregular wave surface receiving kinetic 
energy frcm the airflow can influence the turbulent regime within the 
adjacent airflow. Therefore, empirical formulae, valid or at l eas t 
tested over land, may not be valid over water. 

Examples of the conflicting results obtained frcm strictly empirical 
approaches appear in Figure 1 in which various proposed drag coefficient 
relations are presented for over the water. Being able to add seme 
dynamical description to this region should help in narrowing ones' choices 
for this empirical coefficient. 

Data used in this study were previously analyzed by Davidson (1970) . 
However, those analyses consisted of spectral analyses and the limitations 
imposed on the interpretation of spectral results hindered the description 
of several important properties of tine observed organized motion. Proce- 
dures employed in the present study are well suited for examining organized 
motion within a turbulent regime. The statistical procedures (Holland, 
1968) were designed to reveal systematic structure in various types of 
meteorological data but at the present time have only been applied to 
larger scale fluctuations. This study is the first attempt to apply these 
procedures to small scale fluctuations. 
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FIGURE 1. 



Surmari.es of Drag Coefficients Versus Wind Speed 
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II. THEORETICAL CONSIDERATIONS 



A. WIND WAVE-COUPLING 

In a survey of the understanding of the airflow over water, Stewart 
(1967) forewarned that the air-water boundary is not "just another 
boundary layer", but is unique in that the interface is a part of the 
solution for the boundary layer airflow. The interface is free to move 
and, as observed, it takes advantage of this freedom. Stewart, evaluating 
existing wave generation models (Miles and Phillips) , indicated that in 
1967 there were still insufficient overwater measurements to either verify 
or invalidate existing analytical theories of the interaction between a 
turbulent shear flow and the underlying water surface. He reported, how- 
ever, that even though he examined a very large number of overwater 
velocity spectra, in no case were there sharply defined spectral peaks 
corresponding to the wave spectra peaks. Other studies, however, have 
revealed the existence of wave-induced motion in the air. 

Yefimov and Pososhkov (1969) reported: 

"Laboratory experiments indicate that surface waves induce in 
initially quiet air not only periodic velocity fluctuations but also 
stationary flow in the direction of wave propagation. The velocity 
of stationary flow is considerably higher than would follow if we 
assumed it to resemble Stokes wave flow in the water." 

Volkov (1969) reported, based on interpretations of spectral results: 

"... the wave surface has a considerable influence on all of the 
statistical characteristics of turbulence in the atmospheric boundary 
layer - the fluctuation spectra and dispersion, the turbulent fluxes, 
etc. We therefore succeed in finding seme typical features and 
regularities of this influence by the complex method of investigating 
the processes of interaction between the atmospheric boundary layer 
and the agitated sea surface." 
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Davidson (1970), also interpreting spectral results, concluded: 

"... analyses of the boundary layer over a majority of the earth's 
surface cannot be based on the assumption that the surface is a rigid 
boundary. The influence of the waves on the adjacent airflow has been 
found to be too significant to rely on over-land data to describe the 
boundary layer over waves. Reliable empirical relations giving the 
flux of momentum as a function of the wind profile and density gradi- 
ents will have to be formulated from data obtained over the ocean 
surface itself." 

An analytical model describing coupling between a wave surface and 
the overlying air flew was first proposed by Miles (1957) . In this model, 
the wave-induced perturbation is considered in a coordinate system which 
is stationary with respect to the wave (i.e. , surf -rider's coordinate 
system moving at the phase speed c) . Streamlines in such a coordinate 
system would appear as in Figure 2. The Miles' critical height (dashed 
line) occurs at the level where the wind speed is equal to the wave speed. 
The wave, Figure 2, is stationary in this coordinate system; therefore, 
below the critical level the air appears to move to the left since it 
moves slower than the waves. Above the critical level, where the wind 
speed is greater than the phase speed, the air moves to the right. 

The Miles mechanism can be described "as a positive feedback mechanism 
and, as such, it depends on the pre-existence of the wave. Forces within 
the Miles theory leading to wave growth are depicted in Figure 2. 

Region A is a region of negative vortex acceleration and Region B is one 
of positive vortex acceleration, but the absolute value of the vortex 
acceleration is greater at point A than at B. Therefore, there is a 
negative vortex acceleration at that level hence a deceleration of the 
mean wind when averaged over a wave length. This deceleration .results in 
a loss of mcirentum from the airflow which contributes to the momentum 
increase of the wave. 
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FIGURE 2. Streamlines of sinusoidal perturbation in a 
shear flow near Z = Z ; (a) kinematic conse- 
quence, (b) consequence of dynamics of critical 
level . 
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It has been shown (Kinsman, 1965, p. 571) that the kinonatics at the 
critical level, described above, would result in a phase shift, with 
height, of the wave-induced vertical velocity component . 

Therefore, it is observed, within a linear model, that the wave- 
induced fluctuations can be important in altering momentum transfer in 
the air layer adjacent to the waves. Also, specific predictions of this 
theory can be examined, such as a phase shift of the vertical velocity 
with respect to height. It should be noted, however, that within this 
investigation possible non-linear interactions between the wave-induced 
motion and the shear flow will also be examined. 

B. STATISTICAL CONSIDERATIONS AND PROCEDURES 

Statistical procedures used in this study were described by Holland 
(1968). Applied, to three variables, these procedures yield the joint 
probability density function (JDF) for a pair of the variables with a 
conditional mean function (CMF) of a third variable. 

Holland argued convincingly on the suitability of this type of 

analysis for turbulent regimes. In reference to the application of these 

procedures to a turbulent regime, Holland cited Batchelor (1953) : 

"It is a premise of probability theory that a random function f (a) , 
say, defined for all values of a, is determined statistically by the 
complete system of joint-probability distributions of values of the 
function at any n values of a, where n may take any integral value. 
Likewise, the infinite field of turbulent motion is determined 
statistically by the complete system of joint-probability distribu- 
tions of the values of the vector velocity u(x,t) at any n points of 
of space-time." 

A joint probability density function is a means to represent simul- 
taneously observed values of pairs of variables by contours on a two- 
dimensional array. A tri variate statistical relationship is obtained by 
computing a conditional mean of a third variable as a function of the two 
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